(N 

o 

(N 

C 



Oh 

43 



(N 
> 

in 
in 

oo 

en 

o 

(N 



X 
S3 
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We extend our previous work on a derivative expansion for the Casimir energy, to the case of the 
electromagnetic field coupled to two thin, imperfect mirrors. The latter are described by means of 
vacuum polarization tensors localized on the mirrors. We apply the results so obtained to compute 
the first correction to the proximity force approximation to the static Casimir effect. 
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I. INTRODUCTION 

The Casimir force is known to depend on the electromagnetic properties of the relevant objects ('mirrors') and on 
their geometric configuration, in a rather involved way pQ. 

To put the problem we shall deal with in context, let us consider the Casimir force for a quite general situation, 
namely, we assume that the geometry of the problem may be characterized by just two surfaces. Those surfaces 
may correspond, for example, to the boundaries of two mirrors. Alternatively, the surfaces themselves may describe 
zero-width ('thin') mirrors. Yet another possibility is that those surfaces may be the interfaces between media with 
different electromagnetic properties, occupying different spatial regions. In a situation like the ones above, one can 
think of the Casimir energy as a functional of the functions determining the surfaces. Of course, it is generally 
quite difficult to compute that functional for arbitrary surfaces; exact results are available only for highly symmetric 
configurations, the simplest of which is perhaps the case of two flat, infinite, parallel plates. 

However, when the surfaces are gently curved, almost parallel, and close to each other, the proximity force approx- 
imation (PFA) is expected to be a very accurate method to calculate the Casimir energy. Introduced by Derjaguin 
many years ago [2] to compute Van der Waals forces, this approximation consists of replacing both surfaces by a set 
of parallel plates. Then one calculates the energy as the sum of the Casimir energies due to each pair of plates (each 
plate paired with the nearest one in the other mirror). The PFA has also been used successfully applied in other 
contexts, like nuclear physics [3] and electrostatics [4|. 

In spite of the simplicity and long standing usefulness of the PFA, its validity had not been possible to asses until 
quite recently, mostly because there was no systematic way of improving the approximation. Indeed, even the next 
to leading order (NTLO) correction was unknown. In a recent work [SJ, we have shown that the PFA can be thought 
of as an expansion of the Casimir energy in derivatives of the functions that describe the shapes of the surfaces. The 
leading order in this expansion, that contains no derivatives, does reproduce the PFA, while the higher order terms 
contain the corrections. In Ref. [5] we considered the case of a flat surface in front of a gently curved one, the latter 
described by a function X3 = ip(xx, X2). For simplicity, we computed the vacuum energy for a massless quantum scalar 
field satisfying Dirichlet boundary conditions, the result being: 

^"-flio/^!' ^^+^ d ^) 2 ] . W 

with fix = 1 and /3 2 = 2/3. The first term in this expression is the PFA, while the second term is the NTLO correction. 
This result has been generalized by Bimonte et al to the case of two curved, perfectly conducting surfaces, for scalar 
fields satisfying Dirichlet or Neumann boundary conditions, and also to the electromagnetic case [6]. The results for 
the latter are, #1 = 2 and /? 2 = 4/3(1 - 15/tt 2 ). 

As a validity check, it has been shown that, whenever analytic results are available for particular geometries, the 
corresponding derivative expansion correctly reproduces both the PFA and its NTLO correction [SJH]. Moreover, 
initial discrepancies |6] between the improved PFA and the analytic calculations for the particular case of a cylinder 
in front of a plane [7] , has been resolved in favor of the improved PFA after a revision of the rather involved analytic 
calculation for this particular geometry [S]. 

Bimonte et al also considered the case in which the surfaces are interfaces between different media, with frequency- 
dependent permittivity [9]. In this case, the numerical coefficients j3x and 02 become rather complicated functions of 
ij) and of the dimensional constants that describe the electromagnetic properties of the media. 

In this paper, we will extend the improved proximity force approximation to the case of two imperfect thin mirrors. 
This kind of configuration have already been considered in several previous works; for instance, in order to describe 
the interactions of plasma sheets, graphene sheets or, more generally, arbitrary semi-transparent mirrors, both for the 
static and dynamical Casimir effects [TPl - fl"2"] . 

In some derivations of the Casimir energy for perfect and imperfect mirrors, the boundary conditions at the interfaces 
are represented in terms of auxiliary scalar fields coupled to the TE and TM modes of the electromagnetic field [13] . 
We will follow here a similar approach, but developing a new formalism, based on vector auxiliary fields that couple 
to the dual of the Maxwell tensor F a p evaluated on the surfaces. In this formalism gauge invariance is more apparent 
at the different stages of the calculation. Moreover, the formalism could be useful to address problems with more 
complex geometries, where it could not be possible to describe the electromagnetic field in terms of independent TE 
and TM modes. 

The paper is organized as follows. In the next Section we describe the model and introduce the necessary definitions 
and conventions. In Section 3 we derive the formal expression of the vacuum energy for the electromagnetic field, using 
the above mentioned formalism based on a vector auxiliary field. The derivative expansion for the electromagnetic 
vacuum energy is presented in Section 4. We discuss the results in Section 5, where we analyze the two limiting cases 
of perfectly conducting and near transparent mirrors. In the latter, we find that the NTLO correction to the PFA is 



tantamount to use the area of the curved surface in the leading order expression. We also discuss in that section a 
particular class of imperfect mirrors, in which the transmission and reflection coefficients do not involve dimensionful 
constants. For these 'graphene-like' mirrors, dimensional analysis implies that the vacuum energy is of the form given 
in Eq.pJ), where fa and fa are constants and depend on the dimensionless quantities that describe the mirrors. We 
present some numerical evaluations of these coefficients that interpolate between almost transparent and perfectly 
conducting mirrors. Section 6 contains the conclusions of this work. 

II. THE MODEL: DEFINITIONS AND CONVENTIONS 

We shall consider a model in which the role of the fluctuating vacuum field is played by an Abelian gauge field, 
Aft, in 3 + 1 dimensions, coupled to two imperfect mirrors, L and R. These are presumed to have negligible widths, 
so that we shall use an idealized description whereby they are treated as mathematical surfaces. We will, moreover, 
assume that one of the surfaces (L) is a plane, while the other (R), which may be curved, can always be described by 
a Monge patch. Summarizing, the two surfaces correspond to: 

L) x 3 = R) x 3 = ip(xi,X 2 ) , (2) 

where Xi (i — 1,2,3) are the spatial coordinates. Throughout this work we shall use Euclidean conventions, with 
x^ = x^ (p = 0, 1, 2, 3), xq being the imaginary time. However, we have found it simpler to keep our treatment quite 
general regarding the actual form of the surfaces, postponing the use of ^ to the point when we actually need those 
particular expressions. 

The action S for this model will have the following structure: 

S = S(A;y L ,y R ) 

= S {A) + S L (A; y L ) + S R {A; y R ) , (3) 

where A denotes the electromagnetic field and Sq its free action. Sl and S R are terms that couple A to each 
mirror, with yi, and y R denoting parametrizations of their respective surfaces. Sl and S R can be different because 
of two reasons: first, they correspond to different surfaces, and second, they may also have to account for different 
electromagnetic (response) properties, for example, when the mirrors are composed of different materials. 

We shall then consider a rather general term, Ss, corresponding to the coupling to an arbitrary surface E, par- 
ticularizing to the L and R cases afterwards. Thus, we assume the static surface E to be defined in parametric 
form: 

E) (a 1 , a 2 ) -+ y{a\a 2 ) e M. 3 . (4) 

Although the surface is static, to write the (2 + 1-dimensional) term Ss it is, however, convenient, to introduce 
a parametrization for the world-volume V swept by the surface E, since that is the spacetime region V where the 
interaction takes place: 

V) (a , a 1 , a 2 ) -> y» (ct° , a 1 , a 2 ) = y"(er), M = 0, 1, 2, 3, (5) 

where y° = <j°, and y as in Q. The world- volume is three-dimensional, and we adopt the convention of using indices 
from the beginning of the Greek alphabet to denote components in that space; for example, in an expression like da a 
we implicitly assume that a runs from to 2. We do need to introduce more objects in that space, like the induced 
metric, g a p{a), which may be written in terms of the parametrization: 

9«P(°)- daa dQ& ■ (6) 

We also need to introduce e^, a local basis of tangent vectors to V, such that e^ = q^- They are, by construction, 
normalized to satisfy the condition e^^e^Ja) — g a p(a). 

Before writing the explicit expression for the action, and to make contact with previous works, let us describe a 
simpler model with a vacuum scalar field. The free action is 

So(<p) = / d^x-d^ipd^tp. (7) 

Assuming that the surface action Ss is quadratic in ip, its general form is 

Sv{^y) = \ /rfVv^rfV'v^M^^faX^'), (8) 



where g(a) is the determinant of the induced metric and ir(a, a') describes the (nonlocal) response of the mirror. The 
local approximation of this action is 

S s (^y) = ^Jd 3 a^-Ma) 2 , (9) 

where \ s is a constant (the subindex s stands for scalar). Eq.Q describes the so called '^-potentials', widely used 
as toy models to describe imperfect mirrors. One can check that this kind of potentials induce a discontinuity in the 
normal derivative of the scalar field across the surface, i.e. disc[<9„(p] = X s ip. The factor y/g(tj) is crucial to produce 
such boundary condition |14j . In the limit X s — > oo, the field must vanish on the surface in order to have a finite 
discontinuity across the surface, and therefore one recovers the usual Dirichlet boundary condition on E. 
In the electromagnetic case, the explicit form of Sq(A) will be 

So(A) = Jd^F^F^ + b -(d,A^} , (10) 

where the term proportional to b provides the gauge-fixing. In the calculations presented in the next sections, we 
shall adopt the Feynman (b = 1) gauge. 

Assuming that <Ss is quadratic in Ay,, gauge invariance implies that it will have the general form 

Sx(A;y) = ] f d 3 a ^g(ajd 3 a' y/gja^) F af3 (a)F aip ,(a')TT a ^^' (a,a r ) , (11) 



4 . 
where ir a P a ° is a polarization tensor that depends on the microscopic degrees of freedom on the mirror, and 

F aP = V a A p - VpA a = d a A p - d p A a , (12) 

where V Q denotes the covariant derivative operator, corresponding to the connection for the induced metric, acting 
(in this case) on a covariant vector. We have used A a (a) as a shorthand for the components of the gauge field Ay(x) 
on V, projected along the directions defined by the local basis: 

A a {a) = A,[y(a)]e^a) . (13) 

As in the scalar case, we will start our discussion with a local interaction 

Sv(A;y) = - A J ' d 3 a^gJa)F a0 F a ^ (14) 

where A is a constant, and afterwards we shall extend the results to include frequency-dependent couplings. Note 
that the constants A and A s have different dimensions. 

As it should be evident from the actual form of the interaction term that we are assuming for the model, this kind 
of mirror involves only on the gauge field components which are parallel to the world-volume. However, a term like 
this induces discontinuities across the surface of the component of the electric field which is normal to E, and of 
the components of the magnetic field which are parallel to that surface. The discontinuity is proportional to A and 
depends on parallel components of the gauge field, producing the boundary conditions of a perfect conductor in the 
limit A — > oo, as in the scalar case. For instance, for the flat surface at x% = the boundary conditions read 

disc(F 3 ,) = A9 A F^, (15) 

where the sum over fi excludes /i = 3. These boundary conditions can be explicitly written in terms of the field 
components as 

disc(E 3 ) = XdiE x , 

disc(B 2 ) = -Ado-Ei, 

disc(B 1 ) = -A(a ^ 2 + aiS 3 ), 

where, for simplicity, we assumed that the fields do not depend on the coordinate x^. 



Note that in Eq.(14| we are assuming a Lorentz-invariant interaction, which would be produced by relativistic 
degrees of freedom on the mirror. One could of course consider the interaction between the gauge fields and non- 
relativistic matter on the mirror, giving different boundary conditions |15U16j . For example, the boundary conditions 
obtained in the case of a fluid of nonrelativistic electrons [15] coincide with ours for the normal component of the 



electric field, but differ for the parallel components of the magnetic field. Accordingly, the reflection and transmission 
coefficients in both models will be different. 



Finally, one should check that a term like ( f 4 ) does preserve gauge invariance. This is indeed the case that can be 
seen from the fact that under the transformation: A^ix) — > A^x) +d fJ- uj(x), which is a U(l) gauge transformation 
in 3 + 1 dimensions, one has A a (a) — > A a (a) + 5A a (a), with: 

SA a (a) = d M u\y(a)] e£(a) = d^[y(<r)}d a y»(*) = d a u(a) (16) 

where w(<r) = u[y(a)]. Thus Ss is invariant, since 5F a p(a) = 0. 

III. ELECTROMAGNETIC VACUUM ENERGY: AUXILIARY VECTOR FIELDS 

In the functional approach to the Casimir effect, to obtain the vacuum energy, one usually starts from Z, the 
vacuum transition amplitude, or, equivalently, the zero temperature limit of a finite-temperature partition function. 
For the vacuum field A^, in the presence of the two mirrors, the case at hand, Z may be written as follows: 



Z = JvA^ exp[-S (A)-S L (A;y L )-S R (A;y R )\ . (17) 



To proceed, one should integrate the electromagnetic field. To that end, it is convenient to perform first a transfor- 
mation of the interaction terms, so that A^ only appears linearly, rather than quadratically. This may be done at 
the expense of introducing auxiliary fields, a procedure that we implement now. To simplify the procedure, we first 
represent F a p in terms of its dual F a , a pseudo-vector, such that: 

F a p = VflM^f 1 , \fg~F)F a = e a ^dpA 7 , (18) 

where we adopted the convention that e Q/ g 7 , as well as e a "^ , denote the Levi-Civita permutation symbol (i.e., without 
including any power of g as a factor). 
The generic term Ss may be written as 



Sz(A-y) = ± J ' <P<jJg{aJF a F c 



(19) 



Then we introduce a pseudo- vector auxiliary field £ a (a) , so that the exponential of the interaction term above may 
be obtained as the result of a Gaussian integral: 

exv[-S s (A;y)] = ^- j V$ exp [ - S q {t, A) + ij d 3 <j^)U°)F a (<j)] , (20) 

where S g (£;A) = i ffPay/gWUWi*) and 

AT 6 = fv^e- s ^ x) . (21) 

Note that the representation above is not unique, in the following sense: defining the longitudinal (I) and transverse 
(£) components of £: 

SV) = V^V^' 3 , ff(cr) = £» - £» , (22) 



we see that £/ does not couple to A a . Indeed, because of (18), we see that: 



<P<Ty/g(r)Z a {<T)F a (<T) = d 3 ay/gJ^tf(v)F a (cT), (23) 



where we have used Bianchi's identity: V Q F a = e a,3l d a dpA 1 — 0. 

It is, therefore, possible to modify the auxiliary field action, for example by adding a term depending only on £; to 



the S q term, such that: 



S q -j. S'= S q + % , (24) 



where S^ E is a function of £ (and not of £t) whose precise form will be determined in order to simplify the calculations. 
Thus, a more general (but equally valid) way to rewrite the interaction term is 



m 



exp [ - S^{A- y)]= Tri / n exp - S' q (& X)+i d\^(a)^ a {a)F a {a) 



(25) 



A% = iv£e 



S' q {Z\\) 



(26) 



Besides, note that the term which couples linearly the gauge field to the auxiliary field, can be reinterpreted as an 
interaction with a surface-dependent 'current' j£(x): 



d 3 a^ a (cr)e a ^d A 7 (a) = / d 4 xj£{x) A M (x) , 



where 



J£(a0 = d 3 ad^[x-y(a)]e^a)e^d^(a), 



(27) 



(28) 



is a 'topologically conserved' current, namely, it satisfies 9 M J^ = 0, by its very form, regardless of dynamics. 

The process introduced above for £ may be then independently applied to the two interaction terms, Sl and iSr, 
which are defined as follows: 



S L {A-y L ) = S^(A-y) 
S R (A;y R )=Sj:(A;y) 



S-s-L, y—>-y L , A-s- A L 



E-tR,y-ty R , \-y\ R 



(29) 



Since we introduce one auxiliary field for each interaction term in the action, the gauge field will be coupled linearly 
to the sum of two currents. Indeed, the use of two auxiliary fields, £l and £r, in the partition function, yields: 



A/; 



l"t* I 



VZlDZr e 



S' q (U\^L)-S' q ^ R ;X R ) 



■ DAe -So(A)+ l J d A x J^(x)A^x) 



where J = Jl + J R , with Jj, and J R obtained from (28), replacing £ by L and R, respectively. 
Integrating out A^ in (30), yields: 



3 (a;At)-5' g (5 H ;A H ) 



Z = NTk-J V ^R{e- S '< 

x exp [ - l - J d 4 x J d 4 x'J»{x)D^ (x, x')y (x')] } , 

where D^ is the (free) A-field propagator, which in the Feynman gauge becomes: 

D^/ (#, a/) = Sft/t' D(x, x') . 

Here D(x,x') is the Euclidean free scalar field propagator in 3 + 1 dimensions: 

1 f d 4 k e ik -( x ~ x ') 

D(x,x') = (x\-^\x') = J — — ^, 



(30) 



(31) 



(32) 



(33) 



where we have used a 'bra-kef notation to denote matrix elements of functional operators. Note that Zq — J DA e~ s " 
cannot contribute to the Casimir energy, since it is independent of the coupling to the mirrors. On the other hand, 
the normalization factors J\fi J\fi do not contribute either, albeit for a different reason: each one of them depends 
only on the properties of one the mirrors, being adamant to the coupling of the other. Thus, we define the vacuum 



energy, E vac , in such a way that those contributions, irrelevant to the Casimir interaction energy, are subtracted from 
the very beginning: 

E mc = lim (£) , e- r = Z U ^& , (34) 

where T is the extent of the imaginary time interval. 

We then proceed to the evaluation of T, defined in (34). We note that there still remain in this object contributions 
that correspond to mirrors' self-interactions, depending on only one of the mirrors. They will be neglected, since 
our objective is to calculate the Casimir interaction energy between two mirrors, a physical magnitude to which 
self-interaction energies cannot contribute. 

We deal now with the functional integral expression for T, which in view of the above has the following structure: 

r r = /pa%e- Sr ' & ^» (35) 

where 



d A x J d i x'JE(x)D MM/ (x,x')JE (*') 



Sr(a>60 =S q (S L ;\ L )+S q (t; R ;\ R ) 
1 

+ 2 

+ lJ d4x J ^ J t^) D m' (x, x')J& (*') 
+ \j " d A xJ d A x'J&(x)D^ (x, x')j( [x') 
+ \J d^x J d A x'J^{x)D^ (x, x')j£ (*') , (36) 

which is a quadratic form in the auxiliary fields. In order to perform the integral over the auxiliary fields, we need an 
explicit form for the different terms in St- 

Taking into account (pi), we can find the metric tensors and local tangent vector for each mirror; all of these are 
elements that enter in the terms above. In both cases, the parameters a a are chosen as a a = x a , with a — 0, 1,2. 
We also refer to (x a ) as x\\, reserving the notation xm for {x\,X2)- 

For the L surface, the parametrization is then: 

x \\^Vl{x\\), Vl(x\\) = (*||,0) , (37) 

thus, for L we simply have g a p = S a p, e£ — 5^ for /.i = 0, 1, 2, while e^ = 0. For R, on the other hand: 

x\\ -> Vr(x\\), Vr{x\\) = (a;||, V(x||)) . (38) 

Therefore, 

1 

iScp) = | l + (9i^) 2 d^d 2 ^ I , (39) 

d 2 ^d^ 1 + {d 2 ^f 

which implies: J~g = \J\ + (V?/>) 2 - The tangent vectors, on the other hand, are given by: 

«£(*„) = K + ^(x„)^ = { , r + ^(x„)^ h a^i-U. ( 4 °) 

Then we find that: 

j£{x)D^{x,x')j({x') 

x,x' 

#(*„) (d^S aa '-d a 'd' a )D(x ll ,0;x\ l ,0) e(xj,) (41) 



where we have adopted the notations: d a = d/dx a , d' a = d/dx' a , etc. Besides, we have written the integration 
variables as a subindex of the integral. 

For the analogous term that involves the Jr current instead of J^, the corresponding expression is: 



j£(x)D^,(x,x')J£(x') 



x.x' 



£« (*|| ) { [(d^6 aa ' - d a 'd' a )D( Xll , ^(X||); x[ , V(xj,))] 

"fh^'py'drtixw ) & p , v>(4i ) [a 7 a;^(xn , v(x N ) ; 4 , ^(xj,))] }e (4 ) • (42) 



Finally, 



J£(x)D^,(x,x')J£(x') 
£{x\\) (a^^'-^VjD^ii.O^^Cxf!))^^!) (43) 



and 



a:, a:' 



j£(a:)ZV(a:,a;')j£(z') 



Cf(xn) (^^^'-^'^^(^^(xii);^^)^^^) . (44) 

Defining the matrix kernel T, such that 

Sr = \f £(*||)T°a'(*||, 4)6(4), (45) 

■ /a: li. a: li 

where a,b — L, R, the vacuum energy -B vac may be written as follows: 

1 



-E vac = lim 

T^oo 



TrlnT 
2T 



where the trace affects both continuum and discrete indices. 



(46) 



IV. THE DERIVATIVE EXPANSION IN THE ELECTROMAGNETIC CASE 

As already stressed, the vacuum energy depends nontrivially on the shape of the R surface, i.e. it can be thought as 
a nonlocal functional of ip. When the R surface is gently curved, almost parallel and close to the L-plane, we expect 
this functional to be well approximated by a derivative expansion: 



rf 2 xn [v cS (^) + z(^){d^y + ...\ 

= £<& + £<& + ... (47) 

In order to evaluate the functions V e g and Z, it is enough to consider a class of surfaces of the form ^(x) = a + ?7(x) 
with 77 <C a. Indeed, for these surfaces, and up to quadratic order in 77, the vacuum energy will be of the form 

E vac ~ J d 2 x,| [V cS (a) + V^(a) V + Z(a)(d jV ) 2 + ...], (48) 

and therefore one can obtain V e s and Z from this expression (note that the term linear in r\ vanishes if a is chosen 
to be the mean value of the distances between surfaces) . Therefore, it is sufficient to perform an expansion of T in 
powers of r/, keeping terms with up to two derivatives of r). 
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Denoting by T^ and TW the order-n terms in the respective expansions for T and T, we see that, up to the second 
order, the expansion for the former is given by: 

r = r<°) + r^ + r< 2 > + ... (49) 

where the zeroth and first order terms are: 

r (o) = i Trln [ T (o)] 

r (D = 1- & [(t(°))- 1 t< 1 >] , (50) 

while the second order term receives two contributions I^ 2 ) = r' 2 ' 1 ) + T^ 2 ' 2 \ where: 

r^i) =lTr[(T(°)) _1 T( 2 )] 

r (2,2) = _l Tr [( T (0))-l T (l)( T (0))-l T (l)] _ (51) 

Let us now write the matrices T"', for j — 0, 1, 2: 

(ipCi) tpO') \ 
Tg> < ) ■ (52) 

Those matrices are not completely defined until we adopt a specific form for the action S' , which contains an arbitrary 
part that depends on the longitudinal component of the auxiliary field £ a . In order to render the zeroth-order term 
as simple as possible, it is convenient to add the following terms: 

% = ! / d-Z L {x\\)(x\\\-±=\^)d-S L {rf\) 

2 Jx«,x' 2V-<9 2 



** = IJ X x , v / ^v // ^ v ' ?i?(X||)<X||l 27=P l ' T ' l>v '^ (a; ' |) 



(53) 



where we have used the 'bra-kef notation again, this time for a three dimensional space of coordinates. Besides, the 
derivations arc understood also to act on functions defined on this space of coordinates. 

The T(°) matrix elements, which are invariant under translations along x\u may be Fourier transformed: 

T (0) (^|h4) = f -^e ik - { ^- x '\0¥°\k) , (54) 

and the explicit form of its matrix elements, for the gauge-fixing introduced above, is: 

1_ Ifcl 






ffl«^(*) = (i + ir>» 



m] n J.k) = [f &].(*) = |(<W ~ k 4f)e-^ , (55) 



with \k\ = vfc 2 . 

Regarding the terms of order 1 in 77, it is quite straightforward to see that: 



^lI = Kr = , (56) 



so that, after evaluating the terms that mix L and R, the result may be put in the form: 



d k sWxn-x'u) r -\k\a 

(2tt) 3 



(57) 
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Finally, we consider the second order matrix elements. We shall also discard terms involving more than two 
derivatives of r\. Since the Levi-Civita connection involves at least three derivatives of 77, we replace V by d in the 
gauge fixing term S^ R . Thus, this term will contribute to the second order matrix clement RR only through the factor 
depending on the determinant of the metric. 



(2) 
Besides, we see that T^ = 0, while: 



1 



(2)] 
LR1 



<( x h x '\\) 



T 



(2)] 
RL1 



,{X\\, £||) 

d 3 k 



= iG^fW -<9 Q '^)[?7(x||)] 2 j 



ik-( 



ik-{x\\ —x 



;,) 



\k\e 



-\k\a 



Regarding the RR matrix element, we have four different terms: 



where 



T 



(2) 
RR 



I 



(2,1) 
RR 



■T 



(2,2) 
RR 



T (2,3) T (2,4) 
' R RR "r ^RR 



[^RR > } a J^d 



l rl 



^aa'idjVi^w)) 2 ~ ^m^a'j^r/(X||)Sjr?(xj|) <5 (3) (x|| -x\y 



[ T RR] a J X h X '\) = -MfycuS -9a^«)[»?(X||) ~V(^)f 



d 3 k 

(2tt) 3 



,ik-(x«-x\.) 



\k\ 



[T (2 ' 3) l (xu x'n) - € ia ^e ja ' 'dn(yi ll )d , n(K'u) f d k e ik ' {x ^~ x '\\ ] kpkfi ' 
Y^RR \ a ot'\ x h x \\> _ e e ^V^V^jVy^w) J /2w) 3 ^IH 



[« J L> 



\([dMx\\)] 2 +[dMx'0] 2 



2\k\ 
(2tt) 3 |fc| 



(58) 
(59) 



(60) 



A. Evaluation of T (0) 



We recall that r' ' = |Trln [T( )] where the trace runs over all the indices (Lorentz and indices that label the two 
mirrors). To perform that trace it is convenient to note that 

/ J_ + \M \kl e -\k\a\ 
T(0) ( fc )= (fcJ-A 1 . M K(*) 



2 e 

j_ , 1M 

A L ^ 2 





|fe| 

An ' 2 




\k\ )'P\\W, 



(61) 



where we have introduced the transverse (V±) and longitudinal (V\\) projectors, corresponding to the 3- vector fc, 
namely, [P±] (fc) = 5 aa > 



k a k a , 
k 2 ' 



and [V\\] ,{k) — ° fc2 °' . Since these projectors are orthogonal, 

/ J_ |fe| jfcl -|fc|a \ 
\ 2 C Ah + 2 / 



/^ + M 

WM A % 2 _^ + 



1M 
2 



and 



Trln[T (0) ] = 2 x TL S 



1 x TL Z 



d 3 fc 

(2^ 



lndet 



In det 



j_ , 1M 1M»-I*I« 

A ~ 2 2 

i*L-|fc|a J_ , 1M 

2 e A R " r 2 



J_ 4. 1M 
At ^ 2 









1M 
2 



Discarding a-indcpcndcnt contributions, we see that T^> may then be written as follows: 



r (o) 



1 



L 2 T2 



d 3 k 



2 J (2tt) 3 

L 2 TV cfi (a). 



In 



(W 



(£ + f)(£ + J 3 l ) 



-2|fe|o 



(62) 



(63) 



(64) 
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This result coincides with the vacuum energy corresponding to two imperfect, flat, and parallel mirrors separated by 
a distance a, that we had computed previously ([16]) for the particular case Xl = Xr- Note that, given the boundary 
conditions produced by our relativistic model (see Eqs.(14| and (15)), up to leading order the Casimir energy for the 
electromagnetic field is twice the Casimir energy for the case of a scalar field. This was already shown in Ref. |16| . 
where it was also pointed out that for nonrelativistic matter the contributions of TE and TM modes are not equal, 
in agreement with the fact that the TE and TM reflection coefficients are different in this case [15] . 



B. Evaluation of T (1) 



In the previous subsection we obtained the function V c s- Although the evaluation of the term linear in 77 is not 
necessary for our next purpose of obtaining Z, it is useful as an internal consistency check of the calculations. 
Recalling the expression for T 1 - 1 ', we see that we need the inverse of T*- '. In Fourier space, it is given by: 



[T(°>(fc)]- 



1 



1M 
2 



.lM P -l*|o 



D(k) \-k e -Ma ^ + \M 



V\\(k) 



P±(k) 



(65) 



where: 



^>^ + ^ + M*-'*"> 2 



Then, using the notation: A a f, = [T^ ] 



r« = I r 



±RL 



A LR ) aal (x,x')(T^ L ) a , a (x',x) 



This may be evaluated explicitly by introducing Fourier transforms, the result being: 



r (i) = 



TL 2 



d 3 k |fcp 



(2tt) 3 D(k) 



ma /rf 2 xH(x||) 



From Eqs. (64) and (68) one can easily show that 

p(l) = TL 2 Kff(a) | d 2 X||?;(x|| 



as expected from Eq. ( 48 1 



(66) 



(67) 



(68) 



(69) 



C. Evaluation of T (2) 

We present the evaluation of the two contributions to T^ 2 ' separately. 



1. Contribution of T (2 ' 1) 



Let us consider first T^ 2,1 ). It is convenient to recall the form of the inverse of T^-°\ presented in (65), and of T^ 2 \ 
in equations ( 58 ) and ( 60 ) . 



We first note, by explicit evaluation, that the terms [Tl_r]*- 2 ' and [T^i]*- 2 ) in (58) do not contribute to the force. We 

also see that T R ' R in ( 60b can be ignored, since its tensor structure allows it to mix only with the piece of [T^ - 1 ] 
which is proportional to P« . Since neither object depends on a, the corresponding contribution is irrelevant to the 
calculation of Casimir forces. 
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On the other hand, only the 'transverse', i.e. proportional to Vj_, term in [T^ )] must be retained for the rest of 

the terms. Indeed, it is the only part that can produce an a-dependent contribution for the terms Tljj and "^ RR in 

(60). On the other hand, the tensor structure of "^ RR allows it to mix only with the transverse part of [T^ ^] 

The results due the relevant terms, after extracting the term of second order in derivatives, shall have the following 
form: 



^[AkrTS?] = ^TJ^x^x,,)] 2 , (70) 

with b = 1,2,3. 

Interestingly, T^ yields a local term (no need to perform a derivative expansion): 

c(a)= 1 f d 3 k \k^ 

° l(a> 8 J (2tt)3. (JL 1M)(± + M)2r e 2|fc|a ilZ 1 

(where we subtracted a-independent contributions). We could replace |kii| 2 by 2/3|fc| 2 inside the above integral. 



ForTg^, 



However, in its present form, Eq.(71 1 remains valid even in the case in which X R and \l are functions of fcp. 



d f d* Pl _ ,,, ,^ + ¥ 



c 2 (a) = - lini / y^-r-lp + k\\p\ 



fc^oa|k||| 2 y (2tt) 3 ^ "" D{p) 

i r <?k (i-l¥)N 3 

( M )2 > VV 



8 J (2.)3 ( ^ + ^ )( _^ + ^ )2[e2|fc|a 



r l I |fc| w l | |fc| ■> J 

V \ L T~ 2 Aj^t 2 ' 

where again we subtracted a term independent of a. 

Finally, the contribution due to T RR ' is also local, and the result is 03(a) = —02(a). 

2. Contribution of T (2 ' 2) 

The next point is the evaluation of r' 2:2 '. It can be shown, by using symmetries of the matrix elements appearing 
in the expression, that this contribution reduces to: 

p(2,2) = G+U ( 73 ) 

where: 

G = -\Tr(A LL T^A RR T R i) , U = -^(A £R T«A M T«) . (74) 

After some algebra, we see that the result for G may be written as follows: 

1 f d 3 k 



a = 2Jwr m)] ? m (75) 



1 f d p _2a|p+fc| 



GW = -IJ (2^ e [(p-ip + k)) +p\p + k) 2 \d LL {p)d RR {p + k) (76) 



where: 



*LM = %^,d m <p) s ^. (77) 

Note that, as we are considering static surfaces, 77 (fc) is proportional to <J(fco) an d therefore Gik) = G{ka = 0,k||) 
depends only on kii . Up to second order in derivatives, this term yields a contribution with the same form we had for 
r^ 2,1 \ this time with a coefficient: 

^^o^r^- (78) 
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On the other hand, 



u =\Swf mrm) 



(79) 



where: 



U{k) - 4 J ( 27 r)3 e 



(p • (p + A:)) 2 + p 2 (p + fc) 2 d LR (p)d LR {p + k) 



dLB.(p) 



_ \p\ e 



-\p\a 



2 £>(p) 



Up to second order in derivatives, it produces a coefficient: 

d_ 
k-To (9|k|| 



c 5 (a) = lira Q|1 _ A2 U{k) 



(80) 



(81) 



(82) 



The explicit expressions for C4 and C5 can be obtained by computing the derivatives in Eqs. (78) and (82) 



D. The improved PFA in the electromagnetic case 



Using the results of the previous section, we can finally obtain our main result: the improved PFA for imperfect 
thin mirrors. The zero-order contribution to the vacuum energy is: 



E (o) 



cTXii 



d 3 k 
(2tt) 3 



In 



(f) 2 



-2|*#(x,| 



Ui T 2 Aa h 



2 - 



(83) 



Putting together the results of section IV C we see that the NTLO correction to the PFA reads 



6=1 



(84) 



These results can be immediately generalized to the case in which Xl.r become frequency dependent, for which the 
answer is obtained by making the replacement Xl,r — > X,L,R{ko) i n the final expressions for the derivative expansion 
of the energy. This simple generalization is valid because we are considering static surfaces: the time-dependence 
of the modes of the vacuum field is trivial and therefore one can treat independently each frequency. By the same 
reason, the generalization to the case in which the electromagnetic response of the mirrors depends nontrivially on 
k|| is not so trivial, and would involve local momentum expansions at the different points of the curved surface. This 
point deserves further investigation. 



ANALYSIS OF THE RESULTS 



In this section we perform a general analysis of the results obtained. Let us first consider the perfect-mirror limit. 
One can derive the idealized limit of perfect conductivity by taking Ai, A2 — > 00. The leading term reads, in this case 



E (o) 



pcrf 



(Txii 



d 3 k 



In 1 



-2|A#(xi,)l 



7T~ 

720 



V> 3 



(85) 



which is the well known PFA for perfect conductors. 

The evaluation of the NTLO is more tedious. From the explicit expressions of the functions C\, c-i and C3 presented 
in Section TlV C 11 we obtain: 



c lporf — 
^2pcrf ^3pcrf 



C(3) 
12tt 2 V> 3 



(86) 
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It is worth to remark that, in this limit, Ci per f + C2 pcr f + C3 pcr f = 0, that is, there is no contribution from r*- 2,1 ' for 
perfect conductors. 

On the other hand, the computation of the functions C4 and C5 defined in Section [IVC2| is less straightforward. We 



write the integrals in Eqs. (76 1 and (80) in spherical coordinates in momentum space, perform the derivatives with 
respect to |k||| 2 , and finally compute the integrals. In this way we obtain 

_ 15 - 7T 2 
C4pcrf - C 5 porf ~ J^Q^3 ' ^' > 

and therefore Z(ijj) = j§§&3, in agreement with the result obtained by Bimonte et al [5]. 

It is also of interest to analyze the opposite limit, in which the mirrors are almost transparent ( Al, \r <C ip). The 



zcroth order vacuum energy can be easily obtained by expanding the argument of the logarithm in Eq. ( 83 ) . The 
result is 



The evaluation of the coefficients Cf, is also rather simple in this limit. For instance, from Eq. (71 1 we obtain 



c 1 (V) = i/^AHA i |k || | 2 e- 2 l fe l^ii). (89) 

As already mentioned, c^ + C3 = 0. For almost transparent mirrors, the function D{p) defined in Eq. (66) becomes 
D(p) « (AlA^) -1 . Therefore, up to quadratic order in Xl.r we have 

CiW = -4/(2^ A ^( |fc|2+ 2 |k " |2 ) e "' fe " (X " ) ' (90) 

and C5 = 0. Combining these results we obtain 

E& = ~ljd\ [ djri ( Xl] )] 2 J ^ \ L X R \k\ 2 e->^0 , (91) 

and the improved PFA is therefore 



£vac = -7 / rf 2 x,| V9 I 7^3 X L X R \k\ 2 e-^^0 . (92) 



1 f , 2 r f d 3 k 

Note that, in the semitransparent limit, the NTLO correction is equivalent to the insertion of the factor ^/g, that is, 
the improved result corresponds to the use of the area of the curved surface in the usual PFA. It would be interesting 
to check if this result is exact, as is the case for a scalar vacuum field [T7]- It is also worth to note that the scaling 



of the energy with ip depends of course on the choice of Xl.r'- if both are constants, the integrand in Eq.(92) is 
proportional to ^~ 5 . 

We end up this section with a discussion on the case of graphene-like materials, in which the electromagnetic 
properties of the mirrors are described by dimensionless quantities. As is well known, the charged degrees of freedom 
in graphene can be effectively described by massless fermions confined to the surface, with a propagation velocity 
vp sa 1/300. When the surface is flat, the interaction between charges and the electromagnetic field is described by a 
vacuum polarization tensor which, in our notation, would correspond to a nonlocal interaction with X(k) proportional 
to (fc 2 + k 2 ) -1 / 2 , after a rescaling of the temporal components of all vectors and tensors in order to take into account 
that the fermions have a propagation velocity different from c = 1. For curved surfaces, we expect an effective 
interaction of the form 

5s(A;y) * J d 3 a y/fc) F a ^[V 2 }-^ 2 F^ , (93) 



where V is the Laplacian on the surface (the above mentioned rescaling should also be applied to Eq.(93l). The 
calculation of the Casimir energy for this particular case is of high interest, but is beyond the scope of the present 
paper. We will discuss here a toy model for graphene-like materials, compatible with the assumptions we made so 
far. Therefore, we will consider A£ R = ^L,ij|feo|/2, where £,l,r are dimensionless constants. 
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FIG. 1. Ratio between the coefficient of the leading term Qi(£) and the corresponding value for the perfect mirrors case, oi(0), 
as a function of the dimensionless parameter f . It is a monotonic decreasing function. 
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FIG. 2. Ratio between the coefficient of the NTLO correction 02(C) an d the corresponding value for the perfect case, 0:2(0), as 
a function of the dimensionless £• The inset shows that this ratio is non-monotonous and changes sign for a particular value of 

e- 

By dimensional analysis, the improved PFA will have the same functional form than the case of perfect conductivity, 
that is 



E v 



cTxii 



1 
V>3 



[ai(^L,U) + a 2 (a^R)(d^) 2 



(94) 



The expressions for the dimensionless functions oli can be easily derived from our previous results. We have 
computed numerically these functions for the particular case £l = <tz? = £■ The results are shown in Figs. 1 and 2. As 
expected, these functions approach their perfect conductivity limits for £ — > and vanish in the almost transparent 
limit, for £ — > 00. The leading order in the PFA has a rather simple behavior: the absolute value of the coefficient 
a\ is a decreasing function of £ (see Fig.l). The NTLO correction shows a qualitatively different behavior, since it is 
non-monotonous and even changes sign at a particular value of £ (Fig. 2). Moreover, its absolute value falls faster with 
£, that is, for this kind of materials the NTLO correction quickly loses relevance away from the infinite conductivity 
limit. 



VI. CONCLUSIONS 



In this paper we have computed the Casimir energy for thin and imperfect mirrors using a derivative expansion in 
the shape of the surfaces. The leading term in the expansion reproduces the usual PFA, while the term containing two 
derivatives represent the NTLO correction. These results generalize previous works that involved perfect mirrors 0[6], 
and may be regarded as complementary to those for the interaction between thick mirrors [9]. 

The interaction between the mirrors and the vacuum field has been described by a local effective action, which 
is a novel electromagnetic generalization of the ^-potentials usually considered for scalar fields. We also discussed 
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some nonlocal generalizations, which could be useful to describe the interaction between curved graphene sheets. 
To compute the vacuum energy we used a functional approach, and we used an explicitly (electromagnetic) gauge 
invariant approach, whereby the interaction term has been written in terms of vector auxiliary fields coupled to the 
Maxwell tensor. 

We have presented general expressions for the improved PFA for this model, and checked the particular limits 
corresponding to perfect conductors and almost-transparent mirrors. For the particular case of mirrors described by 
a single dimensionless quantity £, we computed the leading PFA and its NTLO correction as a function of £. We have 
found that the NTLO correction has a non-monotonous dependence on £, and that its absolute value drops quickly 
for imperfect mirrors. 

For the sake of simplicity, we considered a gently curved surface in front of a plane. Moreover, we only considered 
the case in which the curved surface can be described by a single function x% — ip(xi, x-i). The results can be extended 
to the case of two curved surfaces described by functions along the lines of Ref. [5] . The generalization to cases in 
which the surfaces cannot be described in this way, for example the case of an object inside another, is far from 
immediate. 
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